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Abstract 

As the Hubbard energy at half filhng is beheved to reproduce at strong couphng 
(part of) the all loop expansion of the dimensions in the SU (2) sector of the 
planar = 4 SYM, we compute an exact non-perturbative expression for it. For 
this aim, we use the effective and well-known idea in 2D statistical field theory to 
convert the Bethe Ansatz equations into two coupled non-linear integral equations 
(NLIEs). We focus our attention on the highest anomalous dimension for fixed 
bare dimension or length, L, analysing the many advantages of this method for 
extracting exact behaviours varying the length and the 't Hooft coupling, A. 
For instance, we will show that the large L (asymptotic) expansion is exactly 
reproduced by its analogue in the BDS Bethe Ansatz, though the exact expression 
clearly differs from the BDS one (by non-analytic terms). Performing the limits 
on L and A in different orders is also under strict control. Eventually, the precision 
of numerical integration of the NLIEs is as much impressive as in other easier- 
looking theories. 
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1 Prologue 



It is a modern achievement that gauge theories and in particular supersymmetric gauge 
theories hide many reahsations of the algebraic geometry theorisation (cf. [1] just as a 
recent monumental reference on the last discovered parallel and many other features). 

More in specific, the AdS/CFT correspondence [2] should be a general dictionary, 
which would equate - among other physical objects - energies of string states to anoma- 
lous dimensions of local gauge-invariant operators of a dual conformal quantum field 
theory. Proving or even testing this duality in full generality may be a formidable task, 
but the integrability properties of the A/" = 4 super Yang-Mills (SYM) theory have 
proved to be extremely useful to understand how it may work and to which extent. 

The identification [3] of the one-loop dilatation operator of scalar gauge-invariant 
fields with bare dimension L with an 5*0(6) integrable chain with L sites, reducing 
in the SU{2) subspace to the spin 1/2-XXX Heisenberg chain, allowed, by using the 
Bethe Ansatz technique to test the one-loop AdS/CFT duahty in many cases [SIEIE] 
beyond the BMN conditions [H]. In the aforementioned cases, special emphasis has been 
delivered to the 1/L correction, as this would result as the first quantum correction in 
string theory; and more generally all the finite size L corrections would have a similar 
stringy origin and importance. Soon afterwards, integrability of A/" = 4 SYM at higher 
loops started to be hinted and hunted [H]. After various attempts and tests (cf. for 
instance [ID]), eventually in [TT] an all loop asymptotic expression has been proposed 
for the eigenvalues of the dilatation operator in the SU{2) sub-sector, in terms of 
the solutions of Bethe Ansatz-like equations, derived by assuming BMN scaling and 
perturbative integrability. Moreover, this proposal (named after them BDS equations) 
was shown to give the correct (truncated) Bethe equations for the five loop dilatation 
operator, after deriving the latter as an operator. Nevertheless, the BDS equations are 
valid only asymptotically, that is for fixed L when the 't Hooft coupling A is small enough 
that the 0{X^) term (L loops) may be neglected along with the higher order powers. 
In fact, the higher loops are clearly affected by the chain wrapping problem - namely 
an interaction range longer than the chain length - which is not taken into account 
by the BDS proposal. In this respect, an important progress was the remark, by Rej, 
Serban and Staudacher [12], that the SU{2) dilatation operator could be reproduced 
up to three loops by the strong coupling expansion of the Hamiltonian of the half- 
filled Hubbard model. Many tests of the proposal [12] started to be performed (cf. 
for instance [13]), while it seems now clear that starting from four loops the Hubbard 
model will reproduce only part of the entire contributions (likely the 'rational ones'), 
the string theory/gauge theory discrepancies motivating the introduction of a specific 
dressing factor [H] also in gauge theory. Indeed, although the dressing factor would 
also care for the large A behaviour, it is unclear for now how to insert it into the two 
Lieb-Wu Bethe equations for the Hubbard model [T3]. On the contrary, it is manifest 
its introduction into the BDS Bethe Ansatz, and therefore a comparative study of BDS 
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versus Hubbard model is one of the motivations of this paper. 

In a previous paper we proposed a description of the highest and immedi- 
ately lower energy states, for both the 5*0(6) chain and the BDS model, based on the 
non-linear integral equation (NLIE). The NLIE was first introduced in [T7 1 fT8l fT9] for 
studying the finite size scaling of the ground state and of the excited states in (critical 
and off-critical) statistical (lattice) field theories respectively. Although it is equivalent 
to the set of all the Bethe Ansatz equations, it is often more suitable for numerical 
and analytical calculations, especially when it is important, like in the present case, to 
detect how the anomalous dimension (energy) behaves with the length L (especially in 
large L investigations). In fact, it condensates into a single (or only very few) integral 
equation(s) many algebraic equations. In this respect, we will prove here that it is a 
right tool to deal with the two possible ordering of the limits of large size L and large 
coupling A. Furthermore, we will find a systematic way to perform the two expansions 
for small coupling and large coupling at any fixed size. 

In this paper we want to introduce the NLIEs as a profitable treatment of the 
Hubbard model, especially for the understanding of the exact scaling behaviour of the 
dimension (energy) with A and L. We will concentrate on the highest energy (anomalous 
dimension) state of the half-filled Hubbard model {SU{2) sub-sector of A/" = 4 SYM). 
This state is described by two coupled NLIEs which will be written in Section 3. In 
Section 4 we will give an exact expression for its energy, as a function of the coupling 
g and the length of the chain L, in terms of the solution of the two coupled NLIEs. 
This peculiar expression for the energy allows a comparison at large length L (Section 
5) with the analogous result coming from the BDS chain: we will show that the 1/L 
leading term and all the next finite size corrections (power-like and logarithmic) in fact 
coincide (i.e. the usual large L asymptotic expansion do coincide), the difference being 
captured by exponentially small corrections (whose leading contribution we estimated 
at strong coupling). Moreover, as a consistency check of our findings, in Section 6 the 
weak and strong coupling limits of the NLIEs on one hand and of the energy on the 
other hand will be studied and shown to reproduce the known results. The strong 
coupling for the BDS model is also carefully analysed. As a consequence. Section 7 is 
devoted to the understanding of the ordering of the two distinct limits A +oo and 
L +00, both in the Hubbard and BDS models. Eventually, a detailed numerical 
analysis is carried out in the last Section 8. 

2 The Hubbard model: a bird's-eye view 

The Hubbard model was introduced as a simplified model for strongly correlated elec- 
trons on a lattice [20] • In one dimension, it describes A'^e electrons moving on a chain 
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with L sites and interacting via the Hamiltonian 



where cj^, Ci^a are (fermionic) canonical creation and annihilation operators respectively, 
t is the strength of the kinetic nearest-neighbour hopping term, U the coupling constant 
of the density potential and, for our interests, periodic boundary conditions are assumed, 
i.e. Cj_|_x,^(j = Cj^o-) ^i+L.o- ~ ''1,0-- 

In the relevant paper [12], a precise sub-set of the energies E of fl2.1|] was conjectured 
to be proportional to the anomalous contribution 7 to the conformal dimensions in the 
SU(2) scalar sector of A/" = 4 super Yang-Mills theory in the planar limit, 

7=^E. (2.2) 

provided we restrict ourselves to the half-filling case = L and also equate the length 
L to the number of constituent operators. 

Now, a very important part of the correspondence between an integrable system 
and a gauge theory is the mapping of the coupling constants, and the latter can be 
easily argued to be a strong-weak coupling duality for many reasonfl Therefore, one 
possible choice, reproducing the known results up to three loops, may well be 



1 27r 1 Stt^ ^ ^ 

where A = A^fi-yA/ = Sn'^g^ is the 't Hooft coupling of the SU{N) SYM theory in 
the planar limit (A^ — > 00). But this can be modified by higher order contributions 
still preserving, of course, the matching outcomes. Actually, to have a loop expansion 
of fl2.2l) in g"^ with the right wrapping phenomenon occurring at 0{g^^), we need to 
introduce a constant magnetic flux |12j . 



L L 

distinguishing odd and even lengths: = when L is odd and <P = ^ when L is even. 

The Hubbard hamiltonian (12.11) describes an integrable model (infinite many con- 
served charges in involution), which was diagonalised by Lieb and Wu by Bethe Ansatz 
in 1968 [T3]. The twisted Hamiltonian (12. 4p is still integrable and the Lieb-Wu equations 



simple one is the limiting case of the Hubbard model in strong coupling (and half filling) [2T] . 
i.e. the 1/2— XXX revealed, originally, at one loop by Minahan and Zarembo [3]. 
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easily generalise. In the half-filling case they read 



-j-j- M/ - f sm{kj + 



\ui- ^ sm{kj + </)) + 



2 



ui - ^ sm{kj + </)) + ! TT 



^in,-|fsin(fc, + 0)-| " 

where M is the number of down spins. The spectrum of the Hamiltonian is then given 
in terms of the pseudo- momenta kj, by the free dispersion relation 

L 

E = -2t^cosCkj + . (2.6) 
i=i 

Starting from here, we will equivalently derive two coupled nonlinear integral equations 
(NLIEs) for the antiferromagnetic state of the model at any value of L, thanks to the 
methods used in [TB]. For reason of completeness, we point out that the thermodynamics 
(infinite length L, but finite temperature) of the Hubbard model has been studied 
[22l [23] by means of three NLIEs (for a summary of the procedure and a complete list 
of references see [21] )• This approach was based on the equivalence of the (quantum) 
one-dimensional Hubbard model with the (classical) two-dimensional Shastry model. 
For the gauge theory understanding, we need to obtain energies of the Hubbard model 
at zero temperature, but at any value of the length L. This completely justifies our 
approach, and a fortiori in the perspective of extending our calculations to excited 
states (i.e. lower dimension operators in the SYM spectrum). 



3 Two non-linear integral equations (NLIEs) 

Looking at the Bethe equations (12. 5p . we define the function 

<D(x,0 = ^ln^f^, (3.1) 

with the branch cut of ln(2;) along the real negative 2;-axis in such a way that — vr < 
arg2; < tt. Then, we perform a gauge transformation which amounts to adding the 
magnetic flux: 

kj = kj + (l). (3.2) 
After a possible choice of the counting functions as 

W{k) = L(A:-0)-^<|.(^«,-|sinA;,0 , (3.3) 

^ f 2t 1\ 
Z{u) = J]$U--sinfc„-J -5^<l>(n-n„,l) , (3.4) 

7=1 ^ ^ m=l 



5 



we can rewrite the Bethe equations, by taking their logarithm, in the usual form of 
quantisation conditions for the Bethe roots {kj,ui}, 

W{kj) = ■kM + 2-kIJ, (3.5) 
Z{ui) = 7r(Af - L + 1 + 2/f ) . (3.6) 

From now on, we specialise our treatment to the highest energy state, consisting of 
the maximum number M = L/2 of real roots ui and of L real roots kj. For simplicity 
reasons, we restrict ourselves to the case M G 2N (the remaining case M G 2N + 1 is a 
simple modification of this case), which obviously implies L G 4N. 

In the definition of the counting functions (13.31 13. 4p we have to deal with sums 
of functions computed on real Bethe roots, kj and Ui. Let us first concentrate on 
functions of kj. We notice that kj may run only within the first Brillouin zone [— vr, vr) 
and that the functions of kj involved are periodic with period 27r. On the other hand, 
the counting function W{k) is quasi-periodic on that interval and e*^'-'^^ and W'{k) are 
indeed periodic. Using the Cauchy theorem to circulate the interval [— tt, tt) by a small 
displacement e > (this periodic case has been developed in [2S]), we get 



L 



" dk ^ iW\k-ie)(^^^''-'''^ 

where the two complex integrals along — e < ImA; < e at Re/c = ±7r have been neglected 
thanks to the periodicity properties of f{k) and e*^'-'^-*. A thumb rule to understand 
this logarithmic indicator formula goes as follows: since W'{k) > 00, the first integral 
is simply the logarithmic derivative in the following formula, but the second one is not 
because of the non-analyticity of the logarithm^. Nevertheless, the latter can be simply 
manipulated into a logarithmic derivative of an analytic function plus an extra piece: 



+ I' §-J{t: - fc)^ In [1 - e-"'"->] + ^f(k - K)W'(k - it) . 

To make the last term useful, we can compute it along the real axis without any harm, 
because of the periodicity of /(fc) and W'{k:)\ then, after integrating by parts the two 



^We have numerical evidence for that; in Section 5, we analytically prove this statement when 
L — > oo. 

■^We use the approximation 6*^^^+") = g*w('^)g-eW (fe). therefore, we suppose e <C 1. 
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integrals before it, we arrive at 

^J\k + ze) In [1 - e*^('=+*^)] - (3.9) 

- f ^f'ik - ^e) In [l - e^^^^^)] + f ^mw\k) , 
J-7T ^TTZ Zir 



because the boundary terms vanish as a consequence of the periodicity of f{k) and 



Upon integrating by parts the last term, we finally obtain 



+ 



—f'(k)W(k) + Im / —f(k + te) In [l - e^^('=+^^)l + 
■7r2vr IT 

fik)Wik) 



271 



(3.10) 



We will mainly use such formula in the e ^ 0"'' limit, 

L 



" dk 



f'{k)W{k) + 



TT 



f'{k) Imln [1 _ e^'^^'+^o)] 



/(fc)P^(fc) 
27r 



(3.11) 

which reads as (13. lip because of the supposed analyticity of f{k) on the real axis. 

For what concerns a sum of a generic function g{x) for x being any root Ui (which is 
in principle everywhere in the real axis for the ground stat^, we can go along similar 
steps and repeat the original procedure for x G M [IHl [IH] • In this case the boundary 
terms appearing during the computations can be neglected thanks to different applicable 
reasons. One sufficient set of conditions, which apply to the case g = ^, relevant for 
the derivations of the NLIEs for W and Z, turns out to be 

Z'{±oo + iy) = , —€ < y < e ; g{+oo) = —g{—oo) , Z{+oo) = —Z{—oo) , 
Z(+oo ± ie) = — Z(— oo =F ze) , g{+oo ± ie) = —g{—oo ^ ie). (3-12) 



In formulae we can write 

M 

^9{ui) = - 



1=1 



°° dx 
2^ 



°° dx 

7T 



^g'{x)Z{x) + Im / ^g'{x + ie) In [l + e^^^^+^^^j 



(3.13) 



oo 



or, in the e — limit. 



M 



1=1 



dx f°° dx 

—g'{x)Z{x) + J —g'{x)lm\n[l 



(3.14) 



^ As it is well-known after [THl US] , for other states complex roots and holes have to be included. 
^Obviously, e appearing in the equations for Z is different from the homonymous constant related 
to W. 
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Eventually, we have our building blocks in formulae (13.111) and (13.141) . where the r.h.s. is 
written through each counting function respectively. Let us apply them to the definition 
of the W{k), 

W{k) = L{k-(l))+ —<t>' Ix-^smk,-] Z{x)- 

J^oo 27r \ U 2 J 

_ r ^<^' (a;--smk,-]lm\n\l + e'^^^^''^ , (3.15) 

i-oo TT \ U 2J 

and to the definition of Z{u), 

Z{u) = L$ (n, - | sin k, cos k W{k) - 

_ ^y"^$'|^^_^sinA;,0cosA;Imln[l-e*^('=+*o)] - (3.16) 

_ r ^^'^u- y, l)Z{y) + /" ^^'{u - y, 1) Im In \l + e^^^^+^o)] . 
J-oo 2vr y_oo vr 

Inserting in the equation for Z the expression for W coming from (I3.15p . we get 

/ 1\ 2t r dk , f 2t \\ 
Z{u) = L<l> \u,-j +L- J —<^'\u--smk,-jcosk{k-4>)- 

- r ^'^'(^ - 1) ^(^) + r -"^'(^ - 1) 1^1^ [1 + - 

J~oo 2tI y_oo 



I r ("^ _ IsinA:, cosA; Imln [l - e'^^^+^^)] , (3.17) 
U n \ U 2J 



where we used the following cancellation of terms, 



/"^ ( 2t \\ ( 2t \\ 

dk $' ( X - — sin k, -\ cos k $' [y " sin fc, -J = , (3.18) 

which can be easily proven by performing the change of variable k ^ n — k. We now 
write the equation for Z{u) f l3.17p in terms of Fourier transform^, using 

I.(p,0 = -Pf-V"«l^l, (3.20) 



^We define the Fourier transform f{p) of a function /(x) as given by 

f{p) = / dx e-^P^f{x) . (3.19) 
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where P indicates the principal value distribution. We obtain the following expression, 
- e-\P\Z{p) + 2e-\^\Lz{p)-^e--^ I dk e-''-^''"'' cos k Lw{k) = 

^^e"^ [ dke-''-^''''^cosk Lw{k), (3.21) 



U 

where we used the integral definition of the Bessel function Jo (-2), 

/'^ dk 

_^^..sinfc^ (3.22) 

and also the following shorthand notations 

Lw{k) = Imln [l - e^'^C^+^o)] , Lzix) = Imln [l + e*^(^+^°)] . (3.23) 
The terms proportional to Z{p) are now collected and reorganized as 

and, coming back to the 'coordinate' space, we obtain the first of two nonlinear integral 
equations for our counting functions, 

/°° dn T ( ^^'] r°° 

^ sin{pu) ° ^ + 2 / dyG{u-y) Im In [l + e'^(J^+'°)l - 
00 2p cosh f J^^ 

Of 1 

_ fl dk cos k J -^^ ^ Imln h -e'^('=+'°n , (3.24) 

U cosh {TTu-^sink) L J ' ^ ^ 

where G{x) is the same kernel function that appears in the spin 1/2- XXX chain and in 
the BDS Bethe Ansatz (eq. 2.24 of p]), 

GW = £|e-^. (3,25) 

We notice that the first line of the NLIE for Z (KMi coincides with the NLIE (eq. 3.15 
of [16]) for the counting function of the highest energy state of the BDS model. The 
second line of fl3.24p is the genuine contribution of the Hubbard model. 

We finally remark that NLIE fl3.24p can be written in the alternative form 

— gd (ttu -—sink] +2 J dy G{u - y) Im In [l + e*^(^+*°)] + 

^ r dk f^u-— sin k] Im In [l - e*^('=+*°)l , (3.26) 

J^^ Tx dk \ U J 



after introducing the hyperbolic amphtude (the Gudermannian) gd(x): 

gd(x) = / — = zarctane . 

/n cosh t 2 



(3.27) 



On the other hand, starting from (13.151) and inserting in it the equation for Z (I3.24p . 
we obtain the second of our nonhnear integral equations: 



W{k) 



L 



{k-4>) 



dx 



oo 



dp 



sm 



oo V 

1 



2tp 

IT 



sin k 



Imln [l 



cosh sink — nx) 

At n f2t 2t \ 

— J dh G I — sin h — — sin k j cos h Im In [l — e 



1 + e\p\ 

f e*^(^+^°)] - (3.28) 

iW{h+iO)-^ 



The two equations (13. 24^ 13.281) are coupled by integral terms and are completely equiv- 
alent to the Bethe equations for the highest energy state. 



4 The energy or anomalous dimension. 



The eigenvalues of the Hamiltonian (12. 4p on the Bethe states are given by (12.61) . The 
highest eigenvalue can be worked out by using (13.111) . We get: 



E = -2t 



r ^si^kW(k)- T— sin/tlmlnh -e'^('^+*o)l -l1 . (4. 



1) 



We now insert the NLIE for W (I3.28p and observe the cancellation of the first and last 
terms: 



Therefore, we are left with 



E 



" dk 



dk 



sin k (fc — 0) — 1 = . 



dp , f 2tp 



271 



2t < L — sin k — sin —— sin k 



p 



U 



elPl + 1 



" dk 
— sm k 

. 27T 



dx 



Lz{x) 



2t r dk 



u 



sin k dh G 



cosh sin /c — vrx) 
sm n — sm k) 



U 



cosh Lwih) 



— sin k Lw{k) 



We recognize the presence of the Bessel function 



JAz) 



1 

2TTi 



dk sin k e 



iz sin k 



(4.2) 



(4.3) 



(4.4) 
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in the first three terms of the right hand side (in second and third we have used the 
Fourier representations, e.g. fl3.25p for G). We finally obtain that the highest eigenvalue 
of (12. 4p is expressed in terms of the counting functions Z and W as follows, 



E 




dp Jo if) J. if) 



M + 1 



+ 



dx 



dp 



2tx cosh I 



iJi 



2tp 

IT 



(h) cos h 



dp^i^tj,^-^^. 
i e\p\ + 1 



dh 



IT 



Lw{h) sin h 



wi 



+ E' 



W2 



and E- 



w 



E. 



wi 



E. 



W2- 



(4.5) 



The first line of (14. 5p . namely El + Ez, coincides formally with the expression of the 
highest energy of the BDS chain as given in equation (3.24) of [IB]. However, we have 
to remember that for the Hubbard model Z satisfies a NLIE which is different from 
that of the BDS model. On the other hand, the second line, i.e. Ew = Ey/i + Ew2, is 
a completely new contribution. 



5 The large L expansions of Hubbard and BDS en- 
ergies in comparison 

As it was first noticed by [12], in the L = oo limit (thermodynamic limit) the leading 
term of the highest energy Ebbs of the BDS model coincides with the thermodynamical 
expression of the Hubbard model, the first contribution in (14. 5p . Since we can provide 
exact expressions for energies at any length L, we want to extract more information 
about the difference Ebbs ~ E when L is large, but finite. And we are in the position 
to obtain this for any value of the coupling constant g. 

For the highest energy -Ebds, many detailed results were given in [16], where it was 
expressed as (cf. equation 3.24) 



E, 



BDS 




dp MV2gp)MV2gp) 
p\p\ + 1 



PCO PCO ^ipx 



(5.1) 



(with the usual shorthand Lzj^^q{x 
of the NLIE 



Zbds{x) = L I ^ sin pa; 

, 2p cosh 



: Imln[l + e*^BDs{a;+«o)j-j^ terms of the solution 



p- + 2 j dyG{x-y) Imln [l + e^^^^s^s^+^o)] . (5.2) 

9 J -CO 



We use in this section the parametrization (12.30 and we focus our attention on the 
energy formula (14. 5p . For the purposes of this section, it is convenient to restore a finite 
(but small) value for the parameter e > 0, used in the treatment of the function W. 
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As an effect of that, the last term of (13.241) becomes 

1 



/7r 
dk cos{k + ie) 
■n 



In [1 - e'^^''+'''^] . (5.3) 



cosh [ttu — Tx\plg sin(A; + ze)] 
On the other hand, the last term of the NLIE (13.281) for W takes the form 

dhG V2gsm{h + ie) - V2gsmk cos(/i + ze) In [l - 6^^^^+*^)] . (5.4) 
Finally, -Evvi Eyy^ (14.51) are rewritten as 

j-^ ' \J^oo ^ elfl + 1 J L J 

V2 



r dh 

= —2 Im / — cos(/i + ie) 



W2 



Im — sin(/i + le) In [l - e''^^h+^^)] . 

9 J-n 



All these formulae depend on L through the function In [l — e*^('^+*'^)] . Therefore, we 
have to study such a function when L is large. Since e <^ 1 (see Footnote 2), we can 
approximate, at first order. 



ln[l - e^^('=+-)] = ln[l - e^'^W-^^'W] + ^(e^) . 
If we suppose also that 

eiy'(A;)>l V A; G [-TT, tt] , 



(5.5) 



(5.6) 



(this condition will be better stated in few lines), then the factor exp[— e W^'(/c)] becomes 
very small and we are led to the final approximation: 



ln[l - e*^('=+-)] ~ -e'mk)^-ew'ik) ^^ .^^ 

On the other hand, when L — oo we can approximate W{k) by its 'forcing term'. 



W{k) ~ L 
and, consequently, its derivative by 



(5.^ 



W'{k) ~ L 



dp cos ( V2gp sink) — 

-OO 



(5.9) 



The function in the square brackets has a minimum at k = ±7r, which we call uj{g): 

Jo(x) 



uj{g) = 1 



y_oo l + e\p\ Jo 1 



(5.10) 
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Expanding the denominator in power series and integrating term by term we get 

1 



uj{g) = 1 



poo ^ 

2 / rfxJo(x) V( 

•^0 n=l 



1 



oo 
n=l 



-1) 



n+1 



(5.11) 



1 A- — 



The last expression can be seen as a result of an integration in the complex plane 

dz 1 1 



uig) = 1 + 



1 4- ^ 
^ + 232 



(5.12) 



on a curve F (see Figure 6.4 of Takahashi's book [26]), which surrounds the poles on 
the positive real axis of ^.^^^^ , excluding the origin. We can deform the integration 
contour to the curve consisting of the points 6 + iy, with 6 > fixed and \y\ > p > 0, 
and of a semicircle of radius p around the origin; then we let 6 and p go to zero. The 
pole at z = gives a contribution —1 to the integral in the previous formula. The 
integral on the points \y\ < \/2g is zero by disparity of the integrand. On the other 
hand, the integrand computed for y > \/2g equals the integrand in —y, because they 
contain square roots of complex numbers of the same modulus, but lying just above 
(for y > \/2g) or just below (for y < —\plg) the cut. Therefore we are left with 



dy- 



sinh Try 
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(5.13) 



From this expression, it easily follows that uj{g) > 0. Moreover, one can show that 
ci;(0) = 1 and that 

g ^ 00 ^ uj{g) ~ g exp{— 71-^/2 g) . (5-14) 

We conclude that Vf? the derivative W is everywhere greater than a positive constant: 

W'{k) > L uj{g) > 0, V /c G [— 7r,7r]. As a consequence of this, the assumed conditions 

on e and L can be stated as ^ 

- — --<e<l. (5.15) 
L uj{g) 

We now consider the two l^-depending terms of fl4.5p . E\yi and £"^^25 when 05.151) holds. 
We have the following inequalities: 



wi 



< 2 



" dh 



TT 



cos(/i + ie)\ 



i ' e\p\ + 1 



A{V2gp) 



< 2e"'-^'^(^) 



" dh 



IT 



cos{h + ie) 



i e\p\ + 1 



(5.16) 



The integral contained in this last line is finite, as far as e is sufficiently small: it should 
be sinh e < and this condition is always satisfied, as we will show in the following 
Remark 1. Therefore, we conclude that, in the limit L — ^ 00, 



\E- 



Wi 



</i(^?)e-^'^"^^^ 



(5.17) 
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where we have indicated with fi the function of g (but not of L) appearing in fl5.16p . 
The same conclusion for Ew2, 

\Ew,\< f2{9)e~''^^^^\ (5.18) 

can be obtained, in the hmit L — > cxd, by a similar reasoning. 

On the other hand, the same procedure can be applied to the third term of the r.h.s. 
of the NLIE for Z fl3.24p . which we have rewritten for finite e in fl5.3p . This term marks 
the difference between the NLIE for Z in the Hubbard model and for Z'bds in the BDS 
Ansatz and acts a a forcing term in the NLIE for the difference Zbds — Z. One concludes 
that, in the limit L ^ oo, such a term is exponentially small and, consequently, that 

I^BDs(a^) - Z{x)\ < fz{x,g) e"'"^^^^^ , (5.19) 

with an analogous meaning of the function fz{x,g). 

Now, we turn to the expression for the highest energy in the Hubbard model (14.50 
and discuss its relation with the analogous one (15. ip in the BDS context, when L is 
large. We remark that the second term in the r.h.s. of (14.51) is formally identical to the 
second term of (15. ip . the only difference being that in the latter Z is replaced by 2'bds- 
However, the result (I5.19P implies that the difference between these two terms is indeed 
smaller than fEzid) with /^^ a positive function of g. This finding, together 

with (15.171 15.181) . allows us to state that, for all finite values of g, 

L ^ oo \Ebds -E\< fE{g) e'^^^^^) , (5.20) 

i.e. the difference between the highest energies in the Hubbard and in the BDS model 
is exponentially small at large L. Therefore, not only their leading terms coincide, but 
also all the power-like and logarithmic finite size corrections: this is exactly the usual 
asymptotic expansion for large volume in statistical field theory. As a confirmation of 
this statement, we observe that the 1/L correction to the highest energy of the BDS 
model, found in [16] and expressed in terms of the modified Bessel functions Jq, /i by 

V2 hiV27rg) 

Lng Io{V27rg) 6 ' ^ ' ' 

exactly matches the same result for the Hubbard model, obtained with different methods 
by [271 Elj- This has been studied numerically in Fig. [3 

Remark 1. The variable e > introduced in (13.70 has to satisfy the condition e <^ 1 
(see Footnote 2). In any case, an upper bound for e comes from the condition that the 
integration contour of (13. 7p contains no singularities of the functions f{x) appearing in 
the integrand. As far as the NLIE for W is concerned, the function appearing in the 
integrations is $ (13. 3p . Therefore, the singularities come from terms like 

\og('-±{ui-'^smk)\ . (5.22) 
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More precisely, k is a singularity if 

1 2t - 2t ^ 

- ± [-— Im(sinA;)] = 0, - — Re(sin A;) = . (5.23) 

We concentrate on the first equation that takes the form 

I sinh(Im k) I = ^ > ^ . (5.24) 

At I cos(Ke k)\ At 

The upper bound for e, ej/ > e, is given by the smallest value of Im k, namely 

eM = arcsinh — = arcsinh — — . (5.25) 
At 2V2g 

Remark 2. When g <^ 1, we already know that Ebds - E = 0{g'^^) = C(e2^^°f). 
In consequence of that, statement (15.201) is already known to be valid when g 1. 
The results of this section allow to extend the validity of (15.201) - for the highest energy 
state - also to the non-perturbative region. 

Remark 3. On the other hand, when g ^ 1, we can give an explicit expression for 
the estimated difference ( 15.20p . More precisely, we can exactly evaluate E — Ebbs in 
the double limit L ^ oo, g oo. 

When — s> oo, we have e < = ^ 1- Performing the g ^ oo limit of the 
L oo limit of W{k + ie), we get 

W{k + ie) ~ + 26 + arcsin sin k + ie sgn(cos k)] . (5.26) 

The choice e = eM allows this expression to be an expansion in powers of -, exact up 
to terms 0{l/g). In the same limit, the p-integral contained in the formula for E\Yi 
becomes 

dp ^iV2,psin(h+ie) {^9P) _^ 1 sin/i 
i e\p\ + 1 ^J2g \cosh\' ^ ' ' 

It follows that in the double limit L ~* oo, g ^ oo, 

Ew, = -— Im/" — sin/isgn(cos/i) Inh -e'^('^+™^^"'^)-^^(^+^sncosh)i 
g J-^ TT 

Ew, = -^Im T— sin/ilnri-e*^(^+™^"^'^)-^^(i+^s'^^^^ 
9 J-^ ^ 

Therefore, 

+ Ew, = -2— Im — sin h\n\l- ^ (5.28) 

9 J-k/2 

.v^ dh 



2— e~^"^ / — sin /i sin 2L/i 

9 

2^2 

-e ^9 , 



TvLg 
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where we have kept only the leading term proportional to 1/L and we have chosen 

On the other hand, the term (15. 3p . which marks the difference between Z and ^bds? 
in the double limit L ^ oo, g ^ oo becomes 

/TT 
rfA;COsA;5(sinA;) In [l - e*i(fc+arcsinsinfc)-eL(l+sgncosfc)j _ 
-TT 

= - Im[ln(l-e-2^^) -ln(l-e*^^)] =0 (5.29) 
Therefore, in the double limit L oo, g oo we have Z = Zbbs and, consequently, 

2\/2 

E - Ebbs = Ew, + Ew, = . (5.30) 

This behaviour is typical of the "wrapping effects". A similar results in the context of 
string theory was found in |28j . 

Remark 4. For intermediate values of g we can not make a prediction for the 
"velocity" of the exponential damping at large L of |-Ebds ~ E\. However, numerical 
data in Section are consistent with fl5.20p . 



6 Two limiting regimes: strong and weak coupling. 

Conversely to the previous Section, we want now to explore the Hubbard energy (14.50 
in two limiting regimes, ^ — > and j — for any fixed L. They define, respectively, 
the strong and the weak coupling in the Hubbard model and allow for simplifications 
and comparison between our results and analogous ones obtained by other methods. 
These computations are also useful as tests for the NLIEs of W (I3.28P and of Z (I3.24p . 
Besides, we analyse the analogous limit g +oo of the BDS energy for any fixed value 
of L. 



6.1 Strong coupling limit in the Hubbard model, i.e. weak 
coupling in SYM: large y. 

A well known result of the perturbative expansion of the Hubbard Hamiltonian around 
^ = at (strong) half filling shows that the leading term is the Heisenberg 1/2- XXX 
spin chain Hamiltonian [29]. This Section is devoted to derive how our formalism 
consistently reproduces this result and makes natural a linear expansion beyond this 
order. For this aim it is crucial to observe that the NLIE for W becomes redundant for 
very small ^ = 0. Indeed, the NLIE for Z{x) (13.241) easily reduces to 

= Lgd7rx + 2 /" G(x-y) Imln [1 + e^^(^+*°)] +0 ( ^) , (6.1) 
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and hence it precisely agree with the single NLIE for the spin 1 / 2- XXX chain (equation 
(2.25) of [in]) upon forgetting the O (^) terms, namely 

Z{x) = Zxxx(x) + ■ (6-2) 

In the same limit we evaluate the terms entering the rhs of the NLIE for W f l3.28p . The 
integration term on the first line behaves as follows: 

^ dp sin sink) ^ /2tp\ 2tL , T , 1 / 
L / — : II Jo { ^] = — sink dp- r- + ^ 77^ 

J^^ p l + e\p\ \U J U "^l + eM 

= —sink\n2 + l—j . (6.3) 
The term on the second line contains Lz and is (at least) O {jj), since 

dx hl^E^ = . (6.4) 

, cosh TTX 

The two terms just computed are enough to distinguish the leading order of W, 

W{k) = L{k-^) + (J^^ . (6.5) 

This result can be used in the NLIE (13.241) for Z to show that the third term of the 
r.h.s. is O (^^j; because 

dk COS k Im In [l - e^^C^-'^+^o)] = . (6.6) 

-TT 

Since also the first term of the r.h.s. of (13.241) is O (^^j) can correct (16.21) as 

Z{x) = Zxxx(x) + O (^^^ . (6.7) 
Consequently, the term on the second line of (I3.28P is 

oo cosh sin /c — 7ra;J U cosh irx 

(6.8) 

For what concerns the third line of (13.281) . we use (16. 5p to get 

— - / dhG —{sink -sink) cosh lmln[l-e'^^^+'^^] = (6.9) 

J — TT L 

= -^GiO)J\h cos/iImlnfl-e'^^'^-^+'O)] ^ ^ (^) ' 
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The integral in the last line vanishes, as one can see from the following calculation: 



dh cos h In ■ 



1 _ Q-iL{h-4i-iO) 



dh cos{h + (f)) In 



JL{h+iO) 



\ — Q-iL(h-iO) 



^ _ gjL(/i+iO) 

(i/i (cos/icos(/) — sin/isin0) In ^iL(fa^tO) = ~ = • (6.10) 

The summand containing cos /i cos is odd under the change h —h, so its integral is 
zero. The remaining term is also odd under the change h -k — h, thanks to the parity 
of L. Therefore, we conclude that in the limit jj ^ the solution of fl3.28p becomes 



W{k) = L{k - (p) + ^ sin A; In 2 - ^ sin k 



-^^L^,,,(x)sinhvrx^^ 

.oo cosh Tlx 



f/2 



(6.11) 

Curiously enough, we recognize in fl6.1ip the highest energy -Eafx of the ferromagnetic 
spin 1/2-XXX chain (eq. 2.33 of 



W{k) = L{k- 



2t 
U 



sin k Eafx + O 



ft' 



(6.12) 



We now compute the highest energy 04.51) in that limit. Using (16.21) . it is easy to see 
that the first line of (14.51) . El + Ez, is proportional to the highest energy Eafx of the 
ferromagnetic spin 1/2-XXX chain: 



El + Ez 



-2t 



t 



-Eafx + 0[ — 



U 



f/3 / 



(6.13) 



Among the remaining terms, the first one, Ewi, is 2tO therefore, it can be 

neglected, with respect to El + Ez, in the same limit. 

The last term to evaluate, Ew2, deserves more attention. We evaluate it using (16.121) 
and expanding the logarithm where it appears: 



E' 



W2 



+ 



2t 

2t 

2t 
t 



^ dh 1 1 - (^^(^h+io) 

— sin /i — In 

TT 2i 1 - e-*w/(h-*o) 

^ dfl \ ]^ _ gi-L(fe-<^+iO)+j|f sm(/i+iO)£AFX 

— sin h — In ^ — 57 — ^ 

TT 2i I _ Q~iL{h-4>~iO)-i§sm{h-iO)EAF^ 

^ dh ,1, 1 - e^^(''-'^+*o) 

— sin /i — In TTTT — : — ^+ 

vr 2i 1 — e~iL(h-<f>-iO) 



t/2 ) 



:E 



AFX 



" dh 



sin h 



IT 



JL{h~(f>+iO) 
piL{h-<j>+iO) 



^~iL{h—(f)-~iO) 
p—iL{h—<j>—iO) 



o 



t/2 



Now, we now shift the integration variable h ^ h 



E^ 



W2 



JL{h+iO) 



+ 



2t 

t Eafx 
U iL 



dh . ,x 1 , 1 
— sm m + — m ^-77 — — 



" dh 2,, d , i-e^^e^+^o) 
— sm (h + (pj—ln 



+ 



(- 



(6.14) 
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The first integral is very similar to (16.101) and it vanishes in the same way as fl6.10l) 
does. After an integration by parts, the second integral is brought into the form 



-2t 



-2t 



AFX 



U~iL 

t -5'afx 
U iL 



/ — sm[2(/i + 



— (sin2/icos; 



In 



-iL{h-iO) 



COS 2h sin 20) In ■ 



giL{h+iO) 
a—iL{h—iO) 



The second integral is zero by parity; the first one can be shown to vanish by using 
the change of variable h ^ ^ — h, remembering (Section 2) that L G 4N. We conclude 

that, when ^ — > 0, is 2t C j : so, it can be neglected with respect to El + Ez- 
We conclude that in the limit t/U ^ the energy (14. 5 p behaves as 



E 



-2U- 



2Lln2 



dy 



d 



1 



dy cosh ny 



-2t 



jjEafx 



[/2 ; 



(6.15) 



i.e. the Hubbard energy coincides - except for an overall factor - with the same quantity 
of the ferromagnetic spin 1/2-XXX chain (equation 2.33 of [IS])- Besides, with the 
parametrization (12.31) the overall constant in (16.151) —2t^/U = 1. Therefore, in that 
case we have the exact coincidence 



lim E = -Eafx 



(6.16) 



which actually encodes all the non-linearity of this expansion, being the rest, which we 
omit here, just a linear order by order addition to this. 



6.2 Weak coupling limit in the Hubbard model, i.e. strong 
coupling in SYM: small j. 

On the contrary, we now expand the Hubbard energy around y = in a systematic 
way, though we will stop at the first perturbative order. In fact, by means of the usual 
(time independent) perturbation theory the first terms of the expansion were produced 
by [32] • Instead, here we study the NLIE for W (I3.28P by expanding 

W{k; j) = Wo{k) + jWi{k) +0(7)- (6-17) 
The second term of the r.h.s. of (13.281) is expanded as 

. Jojp) r r°°dp . ,/o(p)^ (U 

— sm(psmA;) ~ ^ — sm(psmA;) hoi — 

00 P 1 + J -00 P 2 \ t 

= L arcsin sin /c + o f y j , (6.18) 



19 



since the order 0{U/t) contribution vanishes: 

dp sm{p sin k) Jo{p) = 0, 



_ EL 
~ IT 

The third term on the right hand side is 



(6.19) 



+00 



dx 



5(sin k) 

2t ^ ' 



Lzjx) _ U 

cosh sin A; - vrx] 2t 



+00 



dx 



00 2^ coshf 



Lzix) 



dx Lz,ix) + 0[ — 



(6.20) 



where Zq indicates the order zero of Z in the hmit U/t — > 0. 
For what concerns the fourth term, we obtain 



^ / dhG 



— sm h — — sm k 



cos h Ly/{h) 



/+vr 
dh cos h 5(sin h — sin k) 
-TT 



jj giWo(/i+iO) 



W,{h) + o{j 



since as a distribution 



2t „ /2t 



hm -G{—x]= -6(x) + 0{ — ] . 



(6.21) 



Therefore, at the order zero in U/t the NLIE for W reduces to 

/+00 
-c 

IT 

dh COS h 6{sin h — sink) Lyy^lh) = (6.22) 
L{k — (j)) + L arcsin sin k — sgn(cos k) [Lwo (k) — Lwo {nsgnk — k)] . 



— sin( p sin k) Jo(p) 
2p 



In the last step we used formula (6.693.1) from |30j . 

Now, it is not difficult to show that the solution of fl6.22p is 

Wo{k) = L{k -(f)) + 27rN{k) , 



(6.23) 



where the function N{k) takes only integer values. As far as energy calculations (14.51) 
are concerned, the knowledge of N{k) is not required. 

Let us now focus ourselves on the NLIE for Z fl3.24p . After the change of variable 
p pU/2t in the integrand of the forcing term, we can rewrite it as follows: 



Z(x) 



L 

2t 
U 



+00 



dp pU X Jo ip) 

— sm TT 

2p 2t cosh2g 



it 



+ CXJ 



dk cosk 



cosh [vrx — ^ sin /c] 
20 



dy G{x — y) Im In [l 
Imln[l-e^^('=+*°)] . 



(6.24) 



When J — > the forcing term is clearly 0{U/t). In order to estimate the last term we 
express the inverse of the cosh function in terms of its Fourier transform and we get 



2t 
U 



dk cos A; 

T 

dk cos k 



dp ipl^^-^sink) L 



oo 
oo 
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dp e 



cosh 2 



Lw{k) 



ipi ^x— sin k ) 



-Lw{k) 



, 27r cosh ^ 

dk cos k 5(sin k) Lwo{k) + O ( — 

V t 



-LwM+Lwo{t^) + 0(j 



(6.25) 



as follows from the form (16.231) of Wo{k). This allows us to say that the solution of the 
NLIE for Z in the hmit U/t ^ is 0{U/t). 

Stepping back to the U/t expansion for W, the results on Z allow to say that the 
term (K^ is in fact 0{U'^/t'^). It follows that, at order U/t, the NLIE for W reads 

n- iWo{h+iO) 



Wi{k) 



dh cos h 6{sm h — sin k) Re: 



'l _ giWo(/i+iO) 

whose solution is Wi{k) = 0. Therefore, we can write that 



W{k) = L{k + 2TrN{k) + o 



t 



(6.26) 



(6.27) 



We are now ready to compute the leading term and its first correction for the energy 
(14.51) of the highest energy state in the limit U/t —>■ 0. For what concerns E^, Economou 
and Poulopoulos re-casted it [21] as an asymptotic series in powers of U/t. We write 
only the first three terms of such series: 

'2L_UI^^ 7LC(3)f/2 
8t 



Er 



-2t 



TT m 327r3t2 
We rewrite Ez after the change of variable p — > pU /2t as 

dp sin 



U 

Ez = 2t— dx 
2t 



t^ 



Lzix) 



(6.28) 



(6.29) 
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Since Imln [l + e«^(^+«o)j jg Q[iJ /t) and the p-integral is also 0{U /t), we conclude that 

Ez = 2tO{Uyt^) . (6.30) 

We are left with the contributions coming from the third and the fourth term of (14. 5p . 
E'vi/^ and E^^, which we rearrange, by reintroducing the function G, as follows 



2t 



" dh 



TT 



2t 



— cofihLwih) I dk — sink G 



U 



— (sin h — sin k) 



2t 



— sm n Lw{h,) . 

TT 

(6.31) 
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Using fl6.2ip . we get 



2t 



" dh 



sin h 



TT 



COS h 



" dh 



— cosh LiVo{h)- — + / — sinh Lwo{h) + o i — 



TT 



u 



(6.32) 

Putting together these two integrals, we are left with the following addition to (16.281) : 



+ = 4t lim 

e— >0 



2 dh 



— sin h Imln [l - e*^°('^+'^)] + o ( y 



IT TT 
2 



(6.33) 



After the insertion of (16.231) in (I6.33p . we are left with the computation of 

h Imln [1 _ e*^('^-'^+")] = 







dh 










TT 


vr 




2 




2t 


TT 

r: 








TT 



sin h In 



]^ Q—iL{h—<l))—L€ 



(6.34) 
(6.35) 



Integrating by parts gives 
2Lt ^2 



TT 



d/i COS h 



^iL{h—ip)—Le 



-iL{h—(f))—L£ 



+ 



]^ QiL{h—(p)—Le ]^ Q—iL(h—<f>)—Le 



(6.36) 



In order to perform these integrations, we can see the ratios involved as sums of geo- 
metric series 



2Lt n 



TT 



dh cos h 



^ ^iLn{h-<l,)-Lne + ^ e" 



'iLn{h—<p)—Lne 



_n=l n=l 

Now the integrations can be easily performed, giving 



(6.37) 



2Lt 



vr 



n=l 



-Lne 



2 COS Ln0 2 COS Lncf) 
1 — Ln 1 + Ln 



Off °° 1 

-— Ve-^'^^ cos Ln0- — ^ . (6.38) 

n=l 



Going to the limit e — and rearranging this series we get 

4 cos Lncj) f U 
E^,+E^, = 2t + 

n=l L'' V 

When < </) < 27r the sum of this series is (relation 1.445.6 of [30] ) 



(6.39) 



E\Y^ + E\Y^ — 2t 



2L cos (f - 
2 ^^^^ 



'U 
o I — 



TT 



(6.40) 



Summing fl6.28p with (16.401) . we conclude that in the limit U/t ^ the energy of the 
anti-ferromagnetic state of the twisted Hubbard model behaves as 



E = -2t 



2 cos (f 



UL 



+ o 



sm 



(6.41) 
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When = we get the highest energy of the Hubbard model at small coupling |32j . 



E = -2t 



TT UL fU 

2 cotan— V o \ — 

L 8t \t 



(6.42) 



On the other hand, according to [12], the Hamiltonian of the twisted Hubbard model 
makes contact with the dilatation operator of the SU{2) sector of TV = 4 SYM if 
= 7r/2L. In this case we get 



E = -2t 



1 UL fU 

+ o 



sin^ 8t \t 



(6.43) 



We could go to higher orders, but this result already matches the findings of [33] ob- 
tained by the usual first order perturbation theory. 

6.2.1 On the strong coupling of the BDS Bethe Ansatz 

In [16] we proposed a NLIE for the BDS Bethe Ansatz and we analysed in a careful 
detail the analytic form of the finite size corrections. Because of the asymptotic nature 
of such an Ansatz, we limited our discussion to the case when the limit L — oo is taken 
first, i.e. for finite g. 

However, in the following we will compare the Hubbard model and the BDS Ansatz 
for any range of and hence it is interesting to return on the subject and discuss the 
structure of the finite size corrections of the strong coupling limit of the NLIEs derived 
in [in], i.e. the residual L dependence when the g ^ oo limit is taken first. 

In the previous sections we have shown that the BDS Ansatz NLIE can be formally 
obtained by those of the Hubbard model by simply taking W{k) = tout court. We can 
apply the same reasoning here and derive the strong coupling behaviour of the energy 
for the BDS Ansatz from the computation of the previous section. 

Hence, by neglecting all those contributions coming from the counting function W{k) 
and collecting only those coming from Z{x), i.e. El + Ez, we obtain 

2LV2 _ A + 7-^C(3)v^ + o(- 
Tig Ag"^ 167r^(^^ \g^ 

Obviously, we have reintroduced the parametrization (12.31) for the constants t and U . 
It follows that the first three terms of (16.441) . all coming from provide, up to order 
0{g~^), the exact large g limit of -Erds for any L. 

One can immediately realize that there is a stark difference between the finite length 
corrections of this expression and those obtained in [16]. We will return to this point 
in the next section. 



9^^^E^^, = -—- — ^ ',V,3 VL. (6.44) 



7 Order of limits analysis 

In order to achieve a satisfactory understanding of the behaviour of the anomalous di- 
mension for any value of the coupling constant, it is important to analyse what happens 
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to our equations when the order of the hmits g,L oo is changed. Such an aspect 
can be addressed both in the twisted Hubbard model (with = ti/2L) and in the BDS 
Ansatz, allowing us to compare them explicitly. 

It is important to stress that the NLIEs derived for those models play a crucial role 
in order to have the sub-leading corrections (in g and L) under control. With this a 
piece of information we will be able to infer some interesting properties about the global 
behaviour of the anomalous dimension. 

The Hubbard model. 
L, (yf — > oo: The analysis of Section [5] allows us to immediately write down the following 
expression in the limit L — oo and fixed g 



V2L f+'^ dpJoiV2gp)Ji{V2gp) , IiiV2Trg) 



2 



9 J-oo V e\p\ + 1 Lug IoiV2ng) 6 ^ ^ 

where we have an explicit expression of the (^-dependent coefficient of the L and 1/L 
terms of the L — > oo expansion. A further expansion in g gives 

9 ^ oo = + ^ _ + . . . 7.2 

TT g oL(7 4g^ 

g,L — *• oo: Let us repeat the same calculation, but reversing the order of the limits. 
Upon exploiting one (known) result of Section 16.21 we easily conclude 

^ ^ oo E = - A + • • • (7-3) 

9 sm^ 4^2 

which gives the exact L-dependence of the 1/g term. By expanding in L we have 

^.„n 2Lv^ 17tV2 L 
IT g 2 bLg 4g^ 

The conclusion is that, for the Hubbard model, the limits commute only at leading order 
in both L and g. The disagreement begins with the first sub-leading correction: it is 
interesting to remark that when the order of the limits is exchanged, such a correction 
conserves the same functional form and the only change is in the numerical coefficient 
in front of it. 

The BDS Bethe Ansatz. 
L,g ^ oo: In ref. |[16] we computed explicitly the large L behaviour of the anomalous 

As explained in Section [5], when L ^ oo the highest energies of BDS and Hubbard 
models coincide, up to exponentially small terms. Hence, expanding in g, we have 
again 

fL,g) 2LV2 TX \f2 L 

^BDS =^ + ^ -^ + --- (7.6) 



dimension which turned out to be 



V2L f+'=^ dpJo{V2gp)MV2gp) , ^2 h{V2ng) 
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g,L — i> oo: This case was discussed in Section [B.2.11 The g ^ oo hmit gives 

ILV2 _ L 7LC(3)a 
IT g Ag"^ IQ-K^g 



^ ^ 00 Ebds = - + l\ + • • • (7-7) 



Therefore, we conclude that 



As expected, in the BDS Bethe Ansatz the order of the hmits commutes only at leading 
order in g and L. It is important to point out that the sub-leading corrections differs 
also in their functional form, because the term which behaves as l/{gL) is absent. 

Remarks. 

1. The discussion of this section has many points of contact with that of Section 3 of 
The main difference is related to the treatment of the sub-leading corrections 



in L for E^^J . If one takes the g ^ 00 limit starting from the NLIEs of [TB], one 



^BDS 

immediately realizes that the sub-leading correction used in ref. [M] (given by 
equation (17. 6p ) is not the correct one, because in this limit the structure of the 
counting function changes dramatically giving the structure observed in (17. 8p . In 
particular our analysis of Section 16.2.11 shows that sub- leading corrections in 1/L 
will appear only beyond the order 1/g^. 

2. As expected from the results of Section [5l if we take first the limit L 00, 
the Hubbard and BDS Ansatz behave the same way. This is because, in such a 
limit the charge degrees of freedom (described by the counting function W in the 
language of the NLIEs) are exponentially depressed. 

3. In the BDS Ansatz it was somehow expected that the limits do not commute, in 
particular because of the so-called " wrapping problem" . What is more surprising 
is that the limits do not commute also in the Hubbard model which is believed 
not to be plagued by such a pathology. 

4. As previously pointed out in [M] , the leading strong coupling behaviour at infinite 
length is the same, no matter what the model and order of limits are considered. 

In the following sections about the numerical analysis at fixed L we will use the following 
expressions for the strong coupling expansion 

^ V2 L 

g sm^ 4 g^ 

^BDS = - — + .... 7.9 

'Kg 4^2 
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8 Numerical analysis 



The equations obtained in the previous sections are suitable for numerical evaluations, 
and, as in [16], it is not difficult to solve them by iteration. Our analysis is mainly 
devoted to investigate the difference between the highest energies (or anomalous di- 
mensions) in the Hubbard model and in the BDS Bethe Ansatz. For this reason, we 
will uniquely use the coupling g and we will use (12.31) to express t and U . In particular, 
we observe the value of the ratio 

For BDS, calculations were also performed in [16]. We remember that this model is 
conjectured to work only when -C 1 and up to the order after which the 

wrapping problem is present. In spite of this we will compare the energies of Hubbard 
and BDS predictions out of the strictly perturbative regime. We will begin with the 
analysis of the Konishi operator (L = 4), then we will study the behaviour of the highest 
energy state in the case of chains with an intermediate number of sites. We conclude 
with a numerical study of the difference between the BDS Ansatz and the Hubbard 
model for large L and finite g in order to provide a numerical support to the analytic 
results of Section 



8.1 A test for the NLIEs: the Konishi operator 

The highest energy state with L = 4, corresponding to the Konishi operator, has been 
extensively studied from the perturbative point of view, both in the context of the BDS 
Ansatz and within the Hubbard model. Furthermore, for the latter, an exact implicit 
expression for the anomalous dimension has been found by Minahan [T^ . 

In the present section we will use such known results as a test for the validity of the 
NLIEs derived in this paper. 

Let us begin with the comparison between our NLIE and the perturbative expansions 
calculated in [12] (for our convenience we pushed the computation up to the thirtieth 
order by using the Mathematica routines in Appendix B of [12], see below). 

The result of such a comparison is summarized in Fig. [H Firstly, it is impressing 
the behaviour of the perturbative expansion when compared to the exact result coming 
from the NLIEs: the convergence seems to be quite slow and the perturbative window 
turns out to be very small. We can also observe in Fig. [T] that the ffist line which 
separates is the one corresponding to (18. 2p and the second is (18. 3p . In the small image 
we have magnified the region where the perturbative expansions and the exact curves 
separate. 

The emergence of such a behaviour can be explained by the appearance of rapidly 
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Figure 1. Comparison of the highest Hubbard and BDS energies (anomalous dimensions) for 
a system with L = 4 sites, corresponding to the anomalous dimension of the Konishi operator. 
The curves indicated by "NLIE" are obtained by the non-linear integral equation - (|4.5|) for 
Hubbard and eq. 3.24 of [H] for BDS - and those indicated by "power" are obtained with the 
power expansions ()8.2I and 18. 3p up to the thirtieth order. The value of g where they get far 
away gives an idea of the convergence radius. The small image is a zoom of the surrounded 
area of the largest one. 

growing coefficients in both the Hubbard model and BDS Ansatz cases 
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^Note that the Hubbard coefficients are (multiple of 6) integers and so are the BDS ones 
properly re-scaled. 
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Hubbard, L=4 



NLIE 
exact solution 




Figure 2. Comparison of numerical NLIE data (the same as in Fig. [T]) and the exact albeit 
implicit solution of [13], from weak to strong coupling. The two curves perfectly overlap in 
the common interval g G [0, 1.5]. 



Another interesting check is given by the comparison of our numerical results with 
the exact (implicit) form for the anomalous dimension of the Konishi operator given 
by Minahan in [13]. As shown in Fig. [2], we found a complete agreement within the 
numerical precision of our computation. This allowed us to estimate our relative error 
to be less than 2 ■ IQ-^. 

It is interesting to notice that, even if L is not at all large, the difference between 
the exact curves for Hubbard and BDS remains small as g is increased. This fact seems 
to suggest that the BDS Ansatz can be considered as a good approximation of the 
Hubbard model even when g is large (i.e. non-perturbative) and L is fixed to a small 
value (and not only in the limit of infinite chain). 

From this point of view it would be nice to study what happens in the strong 
coupling regime. Unfortunately such a regime is difficult to reach with the numerical 
integration of our NLIEs, because of a reduced numerical precision for large g. 

However, it is possible to use the strong coupling expansions eq. (17. 9p for the present 
L = 4 case. We obtain 

3.69552 



E 



E- 



BDS 



9 

3.60127 
9 



9' 



9' 



+ 



15) 



This result is interesting for two reasons. Firstly, it shows that already for L = 4 the 
strong coupling prediction for the BDS Ansatz is a good approximation of the corre- 
sponding result obtained in the Hubbard model. Moreover, one can see that the strong 
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coupling expansions smoothly joins the numerical data in Fig. [TJ we will further com- 
ment on this issue in the next section where we will analyse operators of intermediate 
length. 

8.2 Intermediate length operators 

Since our equations are suitable to the study of the energy at any L, we can use them 
to analyze operators of intermediate length. This is the case where the NLIEs can be 
exploited at their best, because in such an intermediate regime it is quite difficult to 
directly use the Bethe Ansatz equations due to the large number of terms. In particular 
we choose to analyze highest energy states with length L = 12 and L = 40. 




0.5 1 1.5 2 

g 

Figure 3. Comparison of Hubbard and BDS energies (anomalous dimensions) for L = 12. 
The two curves are almost indistinguishable in this range of g and start to separate at the 
right border of the plot. 

Let us begin with L = 12. Figure [3] shows that already at such a small value of L, 
and beyond the perturbative region in the coupling g, the curves computed by means of 
both BDS Ansatz and Hubbard model overlap almost completely. This is an indication 
that even if we are dealing with a chain which is far from the thermodynamic limit, the 
predictions of the BDS Ansatz can be considered as a good effective approximation of 
the Hubbard model behaviour. 

Hence, even if the BDS Ansatz is plagued by the wrapping problem, at a quantitative 
level it is able to reproduce all the significant features of the Hubbard model, beginning 
from L ~ 12. 
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L=12 



Hubbard: NLIE 
str. coupl. 
BDS: NLIE 
str. coupl. 




2.1 2.2 2.3 2.4 



4 

g 



Figure 4. The behaviour of the energies E{g) and i^BDS from smaU to strong couphng is 
plotted here for a lattice of 12 sites. The left branches of the curves are the same as in 
Fig. [3] while the right branches are given by the strong coupling expansions (j8.6|) . In the small 
picture there is a zoom of the region where the branches overlap. 



Again, we use eqs. fl7.9p to describe the strong coupling behaviour 

10.8347 3 
E = - + ... 

i^Bos - 1^4 + -- (8.6) 

9 r 

The behaviour from weak to strong coupling for both the Hubbard model and the BDS 
anstaz is plotted in Fig. HI the left branches are the same as in Fig. [3] (numerical solution 
of the NLIEs), while the right branches are given by the equations fl8.6|] (strong coupling 
expansions). As we stated at the end of previous subsection, left and right branches 
smoothly join. 

Let us remark that it was quite unexpected to find the observed good agreement 
between the Hubbard model and the BDS Ansatz predictions for such a small value 
of L. It is also important to stress the crucial role played by the NLIEs in order to 
obtain the exact behaviour of the energy outside the perturbative domain. As shown in 
the study of the anomalous dimension of the Konishi operator, the perturbation theory 
alone is not enough to reach an overlap with the strong coupling expansion. 
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10 I 1 1 1 1 1 1 

0.5 1 1.5 2 2.5 3 

g 



Figure 5. Comparison of Hubbard and BDS energies (anomalous dimensions) for L = 40. 
The two curves are visibly hard to distinguish in this range of g. The largest reached absolute 
difference is 0.0038 namely a relative difference of 0.019%. 

Finally, we compared the Hubbard and BDS anomalous dimensions for L = 40 in 
the range g G [0, 3]: as expected, the agreement between them is further enhanced and 
the two curves can be hardly distinguished, see Fig. [51 

The interesting feature of this case is that we were able to explicitly follow the 
evolution of the relative difference between Hubbard and BDS from weak to strong 
coupling. We observed that the two curves begin to separate at small g, then they 
achieve a maximum in the relative difference, and after that they start to approach 
again. We think that such a pattern is valid at any L, but the reduced numerical 
precision does not allow us to observe it for smaller values of L. 

Remark. The fact that -Ebds > E at weak coupling, but Ebbs < E a,t strong 
coupling suggests that the curves will cross at some intermediate value of g. This is 
consistent with our numerical observation of the approaching of them as g increases. 
Unfortunately, the reduced precision of our data at large g prevented us to observe such 
a crossing explicitly. 

8.3 Large operators at fixed g 

Another interesting analysis concerns the difference between the highest energies of 
Hubbard and BDS as a function of L and at fixed coupling. Our choices were the 
values g = 1.2, 1.6 and g = 2, which lie in an intermediate region for which our 
asymptotic result fl5.30p does not apply. The solution of our NLIEs gives the result 
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Table 1. Table of data concerning the exponential damping I-Ebds — E\ oc e~°'^^^^ as in 
(|8.7|) and in Sectional Two significant digits are available for the data obtained by numerical 
integration of the NLIEs namely for a{g). The inequality (|8.8|) is satisfied with a slowly 
changing or possibly constant ratio. 
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Figure 6. Difference between the highest Hubbard and BDS energies (anomalous dimensions) 
at different sizes of the system and at two fixed values of the coupling, g = 1.2 and g = 2. A 
logarithmic scale is used on the vertical axis. The linear behaviour is clear indication of the 
exponential damping discussed in Section [5] and summarized in (jS.Tp . Related numerical data 
are provided in Table [TJ 



depicted in Table [Hand in Fig. [6] on a log diagram: for large L [L > 300) the behaviour 
is linear, meaning that the difference between the energies decays exponentially as the 
length L is increased. This confirms our analytical findings of Section Consistently 
with that Section, we introduce the numerical rate of decay a{g) and write 

|£;bds - ^1 oc e-'^(^)^ , a{g)>0. (8.7) 

According to our results of Section 5, the inequality 

euj{g)<eMUjig)<a{g) (8.8) 

must hold and in this respect. Table [U suggests that (18.81) is actually correct for the 
values of g chosen. 
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The values of a{g) shown in Table [T] interpolate between the behaviour at small 
g, a{g) = —2\ng, and at large g, a{g) = Consistently, the values of a{g) in 

Table [T] decrease as g increases. However, understanding how a{g) passes from the 
small coupling to the large coupling behaviour on the basis of numerical data seems 
difficult. Indeed, the comparison of the two plots in Fig. [6] shows that the exponential 
behaviour in L of -E — -Ebds is strongly dependent on the actual value of the coupling. 
In addition, it is clear that linearity (i.e. exponential damping) is reached at values of 
L which rapidly increase with g. 




2000 4000 6000 8000 10000 20000 40000 60000 80000 100000 

L L 



Figure 7. The finite size correction 1/L is compared with the theoretical BDS prediction 
(j5.2ip . represented with the dashed straight Hne. As observed in Fig. [6l for higher values of 
g Unearity is reached at higher values of L. 

Finally, we remember, as pointed out before equation (15.211) . that the exponential 
damping oi E — Ebbs forces the equality of logarithmic-like and power-like finite size 
corrections in the Hubbard and BDS models. The most relevant case is the 1/L be- 
haviour, that is explicitly plotted for the Hubbard model in Fig. [7] and compared with 
the analytical prediction (15.211) coming from the BDS Ansatz. 

9 Summary and perspective 

In this paper we have derived the NLIEs describing the highest energy state of the 
half-filled (attractive) Hubbard model with and without a suitable fiux which is re- 
sponsible for a precise contact with the highest possible anomalous dimension (at fixed 
bare dimension L). In particular, according to the important correspondence pioneered 
in [12], we computed the energy /anomalous dimension for this state/operator, and all 
the other conserved quantity eigenvalues may find an exact expression. The dimension 
is of course a function of the 't Hooft coupling A (and of the operator bare dimen- 
sion L), thanks to a duality connexion between this latter and the Hubbard coupling. 
While exact analytical formula for energy/dimension have been extracted only in the 
strong and weak coupling perturbative regime, numerical solutions of the NLIEs and 
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corresponding energy values can be obtained for arbitrary valueq^ of g and L. In this 
respect, we have been able to provide many plots showing clearly the dependence of 
the highest energy (anomalous dimension) on the coupling g. In particular, we have 
concentrated ourselves on the comparison with the highest energy of the (simpler) BDS 
model. We have shown, first analytically then numerically, that at large sizes {L oo), 
they coincide up to corrections exponentially small (i.e. o{L~°°)), i.e. of non-analytic 
form. Of course, the finite size corrections are also important in the string theory since 
they come out as quantum loop corrections. As regards future perspectives, we would 
like to stress that the NLIE approach is a very good and effective method to provide 
the observables dependence on the model parameters and size. One good quality is that 
the NLIE can be written whenever the Bethe Ansatz equations are available, and even 
under more general circumstances, thanks to its equivalence to some functional equa- 
tions. In particular the application of the NLIE to the so-called string Bethe Ansatz 
[T^ and larger field sectors is particularly desirable, in view of tests of the AdS/CFT 
correspondence. For the time being, indeed, in spite of the progress following the dis- 
covery of integrability in both sides of the duality, such tests could be performed only 
for a limited number of cases, because of the technical difficulty in handling the Bethe 
equations for general length and coupling. In the previous work [16] and in the present 
one we have showed that by means of the NLIE framework we are able to overcome 
this difficulty and provide exact analytic scaling with the length L and the coupling A 
(and also numerical evaluation of the conformal dimensions). 
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